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Abstract 

The  investigation  of  the  performance  of  small  single  PEMFC  was  carried  out  by  employing  a  purposely  designed  test  bench  with  complete  control 
of  the  operational  parameters.  A  MEA  preparation  method  was  also  developed  constituting  the  base  for  testing  new  electrocatalytic  materials  or 
improved  electrode  assembling  techniques.  Beyond  determining  the  polarization  curves,  other  tests  have  been  carried  out  like  electrical  resistance 
measurement  by  the  current  interruption  method  and  voltammetric  characterization.  A  semi-empirical  approach  based  on  a  simplified  mathematical 
model  has  been  used  to  fit  the  experimental  polarization  data,  providing  useful  parameters  for  evaluating  the  electrocatalytic  activity  together  with 
other  significant  data  of  the  MEA  performance.  The  analysis  of  polarization  curves  includes  the  determination  of  the  OCV  by  considering  the 
mixed  potential  at  the  oxygen  electrode  by  fitting  the  initial  data  with  an  arbitrary  oxidation  reaction.  This  can  provide  an  estimate  of  hydrogen 
crossover.  Possible  change  of  ohmic  resistance  with  current  intensity  is  discussed.  The  preliminary  results,  demonstrating  the  feasibility  of  the 
method  are  reported. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Since  the  introduction  of  proton  exchange  membrane  fuel  cell 
(PEMFC)  technology,  whose  significant  advantages  as  power 
sources  started  to  be  recognized  in  the  80s,  a  vast  amount  of 
work  has  been  carried  out  to  understand  the  basic  features  and 
to  optimize  the  design  of  such  promising  system.  The  wide  lit¬ 
erature  available  which  has  been  very  well  reviewed  in  excellent 
papers  [1-4],  provides  the  understanding  of  the  complex  phe¬ 
nomena  governing  the  operation  of  PEMFC.  All  the  energetic 
losses  have  been  extensively  investigated  with  particular  atten¬ 
tion  to  the  kinetics  of  the  electrochemical  reactions,  the  ohmic 
losses  in  both  ionic  and  electronic  conductors,  the  mass  trans¬ 
fer  limitations  of  the  gases.  Many  studies  have  been  carried 
out  by  means  of  electrochemical  methods  with  small  electrodes 
in  laboratory  cells,  mainly  to  investigate  the  catalytic  effect  of 
the  electrochemical  reactions.  These  studies  provide  valuable 
knowledge  of  the  basic  phenomena  but  cannot  be  easily  trans¬ 
ferred  to  describe  the  overall  behavior  of  the  fuel  cell  due  to  the 
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complexity  of  the  membrane-electrode-assembly  (MEA)  sys¬ 
tem.  On  the  other  hand,  testing  of  complete  cells  gives  in  a 
simple  way  the  performance  of  the  entire  system,  but  does  not 
easily  consent  the  separate  determination  of  the  various  ener¬ 
getic  losses.  A  remarkable  research  effort  concerning  new  or 
modified  electrocatalysts  has  been  done  in  recent  years  [5-10] 
devoted  to  the  cost  and  catalytic  activity  issues  in  addition  to 
fundamental  research  [11-15]. 

The  need  for  a  common  basis  of  the  two  different  experimen¬ 
tal  approaches  is  mandatory  when  the  object  of  the  investigation 
is  the  efficiency  of  the  catalytic  layer  of  the  real  MEA,  because 
the  presence  of  important  losses  different  from  those  connected 
with  the  electro-catalytic  behavior  can  mask  the  required  infor¬ 
mation  on  this  last  effect.  To  this  end,  the  evaluation  of  the 
different  losses:  activation  overpotentials,  concentration  over¬ 
potentials,  ohmic  losses,  must  be  obtained  from  the  overall 
polarization  curve. 

Starting  from  these  considerations,  our  group  initiated  a 
research  project  to  investigate  possible  improvements  of  the 
electrode  catalytic  layer  in  the  MEA. 

To  this  purpose,  a  complete  test  bench  for  small  single 
PEMFC  or  stacks  up  to  about  100  W  has  been  designed  together 
with  a  test  methodology  aimed  at  the  evaluation  of  the  catalytic 
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layers  of  various  MEAs  by  separately  determining  the  activa¬ 
tion  overpotential  and  the  other  losses.  A  standard  procedure 
has  been  also  set  up  for  the  preparation  of  MEAs  for  the  testing 
of  new  electrocatalytic  materials  or  new  assembling  techniques. 
A  simplified  yet  effective  treatment  for  the  analysis  of  the  exper¬ 
imental  data  has  been  also  developed,  whose  preliminary  results 
are  reported  in  the  present  contribution.  The  semi-empirical 
analysis  of  polarization  curves  presented  here  is  quite  differ¬ 
ent  from  similar  treatment  found  in  the  literature  [20-24].  Some 
of  the  points  which  render  the  present  approach  different  with 
respect  to  the  cited  literature  are  (i)  complete  Butler- Volmer 
equations  for  both  the  cathode  and  the  anode,  instead  of  Tafel 
approximation,  (ii)  the  open  circuit  voltage  is  computed  as  a 
mixed  potential  by  assuming  an  arbitrary  oxidation  current  at 
the  cathode,  (iii)  an  expression  for  ohmic  resistance  as  a  func¬ 
tion  of  current  intensity  is  introduced  to  account  for  phenomena 
which  require  a  more  sophisticated  treatment,  such  as  increase 
of  membrane  resistance  with  current  due  to  partial  de-hydration 
on  the  anodic  side,  or  non-ohmic  effects  due  to  the  porous  3D 
nature  of  the  electrodes. 

The  approach  presented  in  this  paper  has  been  applied  to 
hydrogen  PEMFC,  but  it  can  be  used  for  other  MEA  based 
fuel  cells,  such  as  DMFC  for  which  the  need  to  improve  the 
performance  of  the  catalytic  layer  is  even  more  stringent. 

2.  Experimental 

The  tests  have  been  carried  out  on  commercial  single  cells 
(ElectroChem  Inc.)  of  two  different  size  (1 .4  W,  5  cm2  and  14  W, 
50  cm2).  The  cells  were  assembled  both  with  the  commercial 
MEA  provided  by  the  supplier  (Nafion®  115  as  the  polymeric 
membrane,  with  a  platinum  loading  of  1  mg  cm-2, 20  wt.%  Pt/C 
on  Vulcan  XC-72  catalyst,  carbon  papers  as  gas  diffusion  layers) 
and  with  MEAs  prepared  in  the  laboratory. 

The  MEA  preparation  process,  derived  from  the  literature 
on  the  subject  [3,16-18],  starts  with  painting  the  catalyst  ink 
on  both  sides  of  the  Nafion®  membrane.  The  membranes  were 
pre-treated  by  boiling  for  30  min  at  80°  C  in  7  M  HNO3  solu¬ 
tions,  5%  H2O2,  1 M  H2SO4  and  distilled  water  consecutively. 
The  procedure  assured  the  complete  removal  of  the  organic  and 
metallic  impurities  which  might  be  present  on  the  membrane 
surface.  Prior  to  the  electrode  preparation  the  membranes  were 
dried  and  flattened.  The  catalyst  ink  was  prepared  by  mixing 
supported  platinum  catalyst  (20  wt.%  Pt  supported  on  Vulcan 
XC-72R,  ElectroChem  Inc.)  and  Nafion  (5%  wt  solution,  Elec¬ 
troChem  Inc.)  in  isopropyl  alcohol  (Aldrich)  and  water  1:1.  The 
weight  ratio  between  Pt/C  and  Nafion  was  7:3. 

The  membrane-electrode-assembly  obtained  was  hot-pressed 
between  two  teflonised  carbon  cloths  (ElectroChem  Inc.)  at  a 
pressure  of  50  MPa  for  2  min,  at  1 15  °C. 

The  PEMFC  test  bench  used  for  the  cell  management  has 
been  purposely  designed.  It  presents,  with  respect  to  the  avail¬ 
able  commercial  units,  important  improvements.  One  of  these 
regards  the  possibility  of  controlling  the  amount  of  water  in 
the  gas  humidifying  process  by  means  of  two  separate  pumps 
dosing  the  water  to  be  added,  after  proper  vaporization,  to  the 
hydrogen  and  oxygen  feeds.  Thus,  the  amount  of  water  in  the 


gases  is  independent  of  the  temperature,  being  limited  only  by 
the  saturation  pressure  at  the  operating  temperature.  Moreover, 
the  complete  water  balance  can  be  determined  by  collecting, 
as  liquid  phase,  all  the  water  coming  out  from  the  cell.  All  the 
settings  of  the  test  bench  concerning  gases  flow,  temperature, 
humidity  are  controlled  by  an  automated  system  connected  to  a 
computer,  which  also  allows  the  acquisition  of  the  test  data. 

Commercial  and  laboratory  prepared  MEAs  were  fed  with 
humidified  H2  and  O2  (usually  100%  RH)  gases  at  atmospheric 
pressure;  the  stoichiometry  flow  of  pure  H2  and  O2  was  usu¬ 
ally  2;  the  operating  temperature  of  the  cell  was  set  in  the  range 
70-80  °C.  All  the  cells  to  be  tested  were  initially  subjected  to 
a  conditioning  procedure  to  equilibrate  the  water  content  inside 
the  membrane.  Two  different  procedures  have  been  used,  the  first 
one,  frequently  reported  in  the  literature,  consists  in  holding  the 
cell  at  a  voltage  of  0.6  V  for  at  least  2  h  until  a  steady  value  of  the 
current  intensity  is  obtained;  the  second  conditioning  procedure, 
which  is  the  one  we  used  for  the  tests  presented  in  this  work, 
consists  in  performing  cyclic,  slow  rate  (0.2  mV  s-1),  polariza¬ 
tion  curves  from  the  OCV  to  0.4  V  and  reverse,  until  (typically 
from  10  to  20  cycles)  the  curve  itself  does  not  change  any  longer. 
For  both  procedures  the  cell  is  operated  with  humidified  H2  and 
dry  O2,  with  stoichiometry  flows  of  2. 

The  value  of  the  ohmic  resistance  of  the  cell  was  determined 
by  the  method  based  on  current  interruption  [19],  employ¬ 
ing  the  electrochemical  AMEL  5000  system  which  allows  the 
recording  of  the  cell  voltage  after  current  interruption  with 
a  sampling  time  as  short  as  10  |xs.  The  apparatus  is  a  com¬ 
plete  galvanostat/potentiostat  with  a  built  in  micro-processor  and 
DMA  features  for  the  recording  of  fast  transient.  This  technique 
provides  a  direct  reproducible  measurement  of  the  electrical 
resistance  of  the  membrane  and  electronic  conductors.  The  only 
limitation  in  the  present  equipment  is  the  maximum  output  cur¬ 
rent  of  1  A.  So  we  use  the  apparatus  for  precise  determination 
of  the  ohmic  resistance  at  low  currents,  while  the  polariza¬ 
tion  tests  are  carried  out  by  another  power  supply  operating 
both  at  constant  current  and  voltage.  A  modification  of  the 
power  supply  system  to  perform  resistance  measurements  by 
the  current  interrupt  method  at  higher  currents  is  under  devel¬ 
opment. 

The  polarization  curves  of  the  cells  were  obtained  by  apply¬ 
ing  a  sequence  of  constant  current  steps,  starting  from  zero  and 
increasing  the  current  at  each  step,  after  reaching  the  stationary 
value  of  the  cell  voltage,  which  typically  occurs  in  a  few  min¬ 
utes  (from  5  to  10)  averaging  the  voltage  values  during  the  last 
minute.  In  some  cases  we  have  also  employed  a  slow  potentio- 
dynamic  technique,  by  varying  the  cell  voltage  at  rates  between 
0.2  and  1  mV  s_  1 ,  the  results  of  numerous  experiments  indicated 
no  significant  difference  between  the  two  techniques,  the  second 
one  allows  to  increase  the  number  of  experimental  points  in  the 
first  part  of  the  polarization  curve,  which  is  important  for  the 
analysis  described  in  the  following  section. 

3.  Analysis  of  polarization  curves 

In  order  to  analyze  the  behavior  of  the  MEA,  particularly  to 
obtain  the  information  about  the  catalytic  activity,  a  simplified 
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model  has  been  developed  by  which  the  computed  polarization 
V(I)  curve  is  obtained  and  compared  with  the  experimental  data. 

This  simplified  model  takes  into  account  the  electrochemical 
kinetics  of  hydrogen  and  oxygen  reactions  by  considering  the 
Butler- Volmer  equation  for  both  processes,  employing  a  refer¬ 
ence  value  for  the  exchange  current  density  for  hydrogen  and 
oxygen  (referred  to  1  cm2  of  geometric  MEA  surface)  derived 
from  literature  data  and  the  values  of  the  transfer  coefficients. 
The  limitations  due  to  mass  transfer  of  the  gas  reagents  to  the 
electrodes  are  introduced  in  the  model  through  the  values  of  the 
limiting  currents  for  hydrogen  and  oxygen,  which  in  turn  depend 
on  the  stoichiometric  coefficients  of  the  gases  under  the  test  con¬ 
ditions.  Moreover  the  value  of  the  electrical  resistance  due  to 
ohmic  drop  in  the  membrane  and  in  the  electronic  conductors  is 
considered.  For  the  polarization  curves  analyzed  in  this  paper, 
the  ohmic  resistance  is  computed  starting  from  the  experimental 
value  obtained  by  the  current  interruption  technique. 

The  analysis  of  the  experimental  polarization  curve  is  carried 
out  in  subsequent  steps  by  considering  the  various  contributions 
according  to  the  relationship: 

V(0  =  £rev  -  (*7act,02  +  *?act,H2  +  r) conc,02  +  'leone,  H2  +  #0(1) 

where  activation  and  concentration  overpotentials  are  consid¬ 
ered  in  absolute  value.  The  computation  takes  into  account  the 
formation  of  a  mixed  potential  at  open  circuit,  providing  the 
value  V(0)  <Eiev,  as  discussed  in  the  following. 

3.1.  Activation  overpotential 

To  compute  the  relationship  between  activation  overpotential 
and  current  density  for  both  oxygen  reduction  and  hydrogen 
oxidation,  the  well  known  Butler- Volmer  equation  is  used: 

i  =  i0  |exp  ^aan-^^-exp  ^-acat^^j  (2) 

where  io  is  the  exchange  current  density  and  aan  =  flz, 
«cat  =  (l  ~  P)z  are  the  transfer  coefficient  for  the  anodic  and 
cathodic  branch  of  the  i  (r))  curve  for  the  given  electrochemical 
reaction.  According  to  the  basic  definitions  in  electrochemical 
kinetics,  f>  is  the  so  called  “symmetry  factor”  and  z  the  number 
of  electrons  transferred  in  the  overall  process.  The  other  symbols 
have  the  usual  meaning. 

In  our  analysis,  Eq.  (2)  is  applied  differently  for  oxygen 
reduction  and  hydrogen  oxidation. 

For  the  oxygen  reduction  reaction,  both  2q  and  ffcat  are 
obtained  by  fitting  the  experimental  data  because,  due  to  the 
extremely  low  value  of  the  exchange  current  density,  even  for 
very  low  values  of  current  density,  the  negative  activation  over¬ 
potential  27  is  sufficiently  high  to  render  negligible  the  first  term 
in  Eq.  (2).  Thus,  the  first  part  of  the  polarization  curve,  up  to 
about  5  mAcm-2,  can  easily  be  fitted  on  the  semi-logarithmic 
plot  V/log  (|/|)  by  the  linear  Tafel  equation:  r]  =  a  —  Mog(|/|), 
where  the  two  fitting  parameters  a  and  b,  which  correspond  to 
the  position  and  slope  of  the  Tafel  line,  provide  the  values  for  io 
and  aCat,  respectively.  Thus,  rearranging  Eq.  (2)  by  neglecting 


the  first  exponential  term,  we  obtain: 

ln| z |  =  In  20  -  —j^rb 

(3) 

RT  ,  .  #7 

0  =  - r;lni0 - In  1*1 

«cat#  «catE 

(4) 

r]  =  a-  b\og\i\ 

(5) 

-2.03  RT 

b=  ,  a=  b  logOo) 

Ol  cat  F 

(6) 

This  is  illustrated  in  Fig.  1;  at  first  a 
used  allowing  the  determination  of  a 

tentative  value  acat  =  1  is 
first  approximation  value 

of  io.  Then,  the  best  fitting  for  both  parameters  is  carried  out  in  a 
few  iteration  steps.  It  is  important  to  note  that  in  the  range  where 
the  Tafel  linear  trend  holds,  only  the  activation  overpotential  for 
ORR  is  considered,  the  other  terms  being  negligible,  so  the  b 
and  r'o  values  for  ORR  are  derived  independently  of  the  other 
fitting  parameters,  which  will  be  obtained  in  subsequent  steps. 

Also  for  the  hydrogen  oxidation  reaction,  Eq.  (2)  is  used,  but 
one  fitting  parameter  only,  the  exchange  current  density,  is  intro¬ 
duced,  while  the  values  of  the  transfer  coefficients  are  taken  from 
literature  data.  This  procedure  is  necessary  because  of  the  lack  of 
a  linear  portion  of  the  experimental  data  attributable  to  the  HOR 
alone.  It  must  be  noted  that  due  to  the  relatively  high  value  of  io 
for  hydrogen  oxidation,  compared  to  the  one  for  the  ORR,  the 
above  procedure  appears  adequate  to  the  scope  of  the  analysis. 
In  fact,  the  estimation  of  the  exchange  current  density  for  HOR 
is  obtained  by  fitting  the  experimental  data  in  the  part  of  the 
polarization  curve  where  ohmic  and  concentration  polarization 
are  negligible.  Fig.  2  illustrates  this  point.  Curve  a  on  this  figure 
is  the  computed  curve  for  the  activation  overpotential  of  oxygen 
reduction  as  previously  indicated,  curve  b  is  the  sum  of  activation 
overpotentials  for  oxygen  reduction  and  hydrogen  oxidation,  this 
last  being  the  difference  between  curve  b  and  a,  the  open  circles 
are  the  experimental  points  (curve  c).  The  open  circuit  volt- 


Fig.  1 .  Determination  of  the  exchange  current  density  io  and  cathodic  transfer 
coefficient  acat  for  the  ORR,  by  linear  fitting  of  the  first  part  of  the  polariza¬ 
tion  curve  corresponding  to  Test  2  reported  in  Fig.  7.  Laboratory  prepared 
MEA,  0.5 mg Pt cm-2,  geometric  surface  5cm2,  T=70 °C,  humidified  H2  and 
02  (100%  RH)  at  atmospheric  pressure,  stoichiometry  flow  2. 
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Fig.  2.  Computed  curves  for  the  activation  overpotential  of  the  oxygen 
reduction  and  hydrogen  oxidation  reactions  for  a  laboratory  prepared  MEA, 
0.5  mg  Pt  cm-2,  geometric  surface  5  cm2,  T=70°C,  humidified  H2  and  O2 
(100%  RH)  at  atmospheric  pressure,  stoichiometry  flow  2.  Curve  a  (Eav  — 
7act,o2)  is  the  computed  curve  for  the  activation  overpotential  of  oxygen  reduc¬ 
tion,  curve  b  ( Etev  —  t)act,o2  +  7act,H2)  is  the  sum  of  activation  overpotentials 
for  oxygen  reduction  and  hydrogen  oxidation,  the  open  circles  are  the  experi¬ 
mental  points  (curve  c).  The  open  circuit  voltage  is  shown  on  the  graph  (OCV), 
while  the  reversible  voltage  (Erev  =  1 . 19  V  in  this  test)  is  not  reported  to  improve 
the  readability  of  the  graph.  Concentration  overpotentials  and  ohmic  drop  are 
negligible  on  this  graph. 

age  is  shown  on  the  graph  (OCV),  while  the  reversible  voltage 
(1.19  V  in  this  test)  is  not  reported  to  improve  the  readability 
of  the  graph.  It  can  be  seen  that  curve  a  (Erev  —  rjacifh  )  departs 
from  the  experimental  data  at  about  5  mAcm-2,  while  curve  b 
(Erev  —  rj act,o2  —  bact,H2)  fits  the  data  up  to  about  15  mAcm-2. 
To  obtain  a  good  fitting  of  this  part  of  the  curve,  accounting  for 
the  activation  overpotential  for  hydrogen  oxidation,  a  high  num¬ 
ber  of  experimental  points  is  necessary  at  the  beginning  of  the 
polarization  curve. 

3.2.  Open  circuit  voltage 

It  is  frequently  reported  in  the  literature  regarding  the  model¬ 
ing  of  single  PEMFC  that  the  values  of  open  circuit  voltage 
markedly  differ  from  those  computed  from  the  equilibrium 
(Nemst)  equation.  It  is  commonly  accepted  that  those  dif¬ 
ferences  are  attributable  to  the  oxygen  electrode,  for  which 
equilibrium  conditions  cannot  be  attained  at  open  circuit.  Since 
the  basic  studies  by  Conway,  Bockris,  Hoare,  Damjanovic  and 
others,  for  a  review  of  the  topic  see  Ref.  [25],  the  reasons  for 
this  behavior  were  interpreted  in  terms  of  parasitic  reactions 
responsible  for  a  mixed-potential  conditions.  Among  the  vari¬ 
ous  parasitic  reactions  which  have  been  assumed,  those  which 
received  major  attention  are  Pt  or  electrode  surface  oxidations 
and  hydrogen  peroxide  reactions.  In  addition  to  these  reactions. 


which  occur  independently  of  the  other  electrode  of  the  system, 
an  important  effect  on  the  open  circuit  potential  is  observed  in 
PEM  fuel  cell,  when  possible  crossover  of  hydrogen  through 
the  membrane  is  considered.  In  this  case,  an  oxidation  current 
for  hydrogen  at  the  cathode  has  been  considered.  The  effects  of 
hydrogen  crossover  have  been  investigated  by  various  authors 
[26-31], 

In  our  simplified  approach  the  above  phenomena  are  taken 
into  consideration  by  assuming  an  arbitrary  oxidation  current 
which  should  include  all  possible  causes.  This  is  schematically 
illustrated  in  Fig.  3,  where  the  fitting  of  experimental  data  at 
low  current  densities  is  obtained  by  using  such  an  empirical 
approach. 

On  this  graph,  the  circles  are  the  experimental  point.  The  fuel 
cell  reversible  voltage  and  the  open  circuit  voltage  due  to  the 
mixed  potential  are  shown.  The  “parasitic  oxidation  reaction” 
curve  is  arbitrary  and  corresponds  to  best-fitting  values  of  the 
equation: 

r,  =  Ho  +  ^i°g(5-y  (7) 

where  Vs  and  /s  are  the  voltage/current  values  for  the  para¬ 
sitic  reaction(s),  Vso  the  half-wave  potential  of  the  curve,  cs  an 
empirical  coefficient  accounting  for  the  curve  slope  and  7is  is 
the  limiting  value  for  the  oxidation  reaction.  Even  if  the  fitting 
by  using  three  parameters  appears  awkward,  it  can  be  observed 
that  the  sensitivity  of  those  parameters  is  such  that  best  fitting 
conditions  can  be  obtained  easily.  Though  this  fitting,  aiming 
at  finding  the  value  of  the  open  circuit  voltage,  is  practically 
independent  of  the  fitting  parameters  corresponding  to  the  other 
dissipation  sources  (oxygen  and  hydrogen  activation  overpoten¬ 
tials,  ohmic  contributions  and  concentration  overpotentials),  it 
is  more  easily  accomplished  if  the  activation  overpotential  for 
ORR  has  already  been  carried  out. 


Current  /  mA 

Fig.  3.  Mixed  potential  analysis  of  the  polarization  curve  at  very  low  current 
density  to  obtain  a  curve  fitting  starting  from  the  open  circuit  voltage.  The 
computation  provides  an  estimation  of  hydrogen  crossover  from  the  value  of 
the  current  of  the  parasitic  reaction  at  the  OCV.  The  data  refer  to  the  same  test 
shown  in  Figs.  1  and  2. 
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An  interesting  result  from  this  computation  is  the  value  of 
the  current  density  for  the  parasitic  reactions  at  the  open  circuit 
potential.  This  current  is  mainly  due  to  the  hydrogen  crossover 
oxidation  and  in  many  experiments  [26,27]  carried  out  with 
Nafion®  membranes  at  70  °C  ranges  from  0.2  to  0.5  mA cm-2. 
These  data  are  in  good  agreement  with  our  computed  values. 

3.3.  Concentration  overpotential 

Concentration  overpotential  arises  from  the  concentration 
gradient  which  is  established  when  a  given  current  flows  within 
an  electrochemical  cell  due  to  the  consumption  of  reagents  at 
the  electrode  surface  and  the  consequent  mass  transfer  (diffu- 
sional)  mechanism  to  replace  them.  The  simplest  way  to  express 
this  type  of  overpotential  is  to  consider  the  value  of  the  limiting 
current  density  corresponding  to  the  flux  of  reagent  supplied  to 
the  system. 

In  case  of  gaseous  reagents  such  as  oxygen  or  hydrogen,  if  we 
know  the  flux  (e.g.  <yn2  mol  s-1  cm-2),  then  the  current  corre¬ 
sponding  to  such  a  flux,  |i/|  =  2  /y/h2  A  cm-2,  is  the  maximum 
possible  value,  i.e.  the  limiting  current  density.  If  we  consider 
concentration  overpotential  alone,  by  neglecting  all  other  dissi¬ 
pation  terms,  the  relationship  between  overpotential  and  current 
density  is  easily  obtained: 


By  considering  the  absolute  value  of  the  overpotential,  such 
a  relationship  is  valid  for  both  cathodic  and  anodic  reactions, 
because  the  logarithmic  expression  is  always  negative. 

If  we  apply  Eq.  (8)  to  the  hydrogen  oxidation  reaction  we 
obtain: 

l^conclfb  =  —  (  1  ~  '.  )  (9) 

2  F  \  H,  H2 ) 

and  for  oxygen  reduction: 

l7concl02  =  — f  1  —  ~  ^  (10) 

4  F  \  z/,o2  / 

This  approach  avoids  complex  calculations  of  concentra¬ 

tion  or  pressure  profiles  of  reagents,  because  it  is  based  on  the 
assumption  that  mass  transfer  limitation  is  described  by  the  con¬ 
centration  gradient  (dC/dx)x=o  at  the  electrode  surface,  whose 
maximum  value  is  proportional  to  the  limiting  current. 

Even  if  the  real  situation  in  3D  electrodes  is  much  more  com¬ 
plex,  the  above  equations  are  applicable  for  the  scope  of  the 
present  work. 

In  the  analysis  of  polarization  curves,  it  must  be  observed  that 
usually  the  flow  of  gaseous  reagents  is  set  to  a  value  described 
by  the  so  called  stoichiometric  coefficient  that  is  the  ratio 
between  the  flow  of  reagent  fed  to  the  system  and  the  flow  cor¬ 
responding  to  the  applied  current.  Since  this  condition  requires  a 
change  of  the  reagent  feed  with  current  to  keep  the  stoichiomet¬ 
ric  coefficient  constant  during  the  polarization  test,  commonly 
the  stoichiometric  coefficient  is  referred  to  the  maximum  value 
of  the  applied  current.  Under  these  conditions,  the  limiting  cur¬ 
rent  value  is:  q  =  A.,imax.  Under  these  simplifying  conditions,  the 


computation  of  the  concentration  overpotential  does  not  require 
the  determination  of  fitting  parameters. 

3.4.  Ohmic  dissipation 

The  computation  of  the  voltage  drop  due  to  the  ohmic 
resistance  of  the  membrane  and  electrodes,  could  be  simply  con¬ 
sidered  as  the  product  RI,  where  the  resistance  value  R,  fitting  the 
experimental  polarization  curve,  fairly  agrees  with  the  measured 
ohmic  resistance. 

It  was  observed  that  in  the  majority  of  polarization  curves 
analyzed,  this  term  is  insufficient  to  provide  a  good  fitting  of 
the  experimental  data,  especially  when  the  test  is  carried  out  at 
rather  high  values  of  current  density  (above  1 A  cm-2). 

It  was  found  that  a  very  good  fitting  can  be  obtained  by 
considering  a  change  of  the  ohmic  resistance  with  current  inten¬ 
sity,  R  =  R(T).  This  assumption  is  consistent  with  the  findings 
of  various  authors  [24,32-34]  who  interpret  the  increase  of  the 
membrane  resistance  with  increasing  current  intensity  as  a  con¬ 
sequence  of  a  partial  de-hydration  at  the  anode  side  due  to 
electroosmotic  drag.  It  is  worth  noting  that  the  change  of  resis¬ 
tance  with  current  is  supported  by  experimental  measurements 
of  the  ohmic  resistance  both  by  the  current  interruption  method 
and  by  electrochemical  impedance  [32,33]. 

We  assume  the  following  relationship:  7?(/)  =  f?0(l  +aln), 
where  Rq  is  adjusted  starting  from  the  experimental  resistance 
value  determined  at  low  current  density  (about  0.2  A  cm-2)  and 
the  parameter  a  is  obtained  by  fitting;  n  is  close  to  1,  but  can  be 
adjusted  (typically  from  0.9  to  1.1)  for  a  better  fitting. 

Fig.  4  illustrates  the  determination  of  the  ohmic  resistance 
by  the  current  interrupt  technique  for  the  laboratory  prepared 
MEA  described  in  the  caption.  The  current  before  interrup¬ 
tion  was  0.998  A,  the  obtained  resistance  value  is  0.24  £2  cm-2, 


T=  70  °C,  humidified  H2  and  02  (100%  RH)  at  atmospheric  pressure,  stoi¬ 
chiometry  flow  2.  The  test  starts  at  a  voltage  of  0.764  V  with  a  stabilized  current 
of  0.998  A.  After  interruption  by  a  fast  solid  state  switch  (<1  p,s)  the  current 
goes  to  zero  and  the  voltage  abruptly  increases.  To  obtain  the  resistance,  the 
voltage  change  is  computed  as  indicated  in  the  figure  by  linear  regression  of  the 
experimental  point  after  interruption.  The  sampling  time  is  10  p.s.  The  measured 
resistance  value  of  0.0485  12  gives  a  specific  resistance  of  0.24 12  cm-2. 
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Table  1 

Results  obtained  by  analyzing  the  polarization  curves  reported  in  Fig.  7 


Polarization 
curve  (Fig.  7) 

Cathode  exchange  currei 
density  (A  cm-2) 

it  Anode  exchange  currer 

density  (A  cm-2) 

it  Cathodic  transfer 

coefficient 

^average 

(Ocm-2) 

Catalytic  activity 
i<t= 0.3  (mAcm-2) 

Current  for  hydrogen 
crossover  (mAcm--) 

Test  1 

8.0  x  10-7 

4.5  x  10-2 

0.94 

0.59 

8.52 

0.18 

Test  2 

5.3  x  10-7 

5.0  x  10-2 

1.0 

0.26 

19.2 

0.28 

Test  3 

1.0  x  10-6 

7.0  x  10-2 

0.99 

0.18 

17.4 

0.34 

which  fairly  well  agrees  with  the  average  value  of  0.26  £2  cm-2, 
obtained  by  fitting  the  polarization  curve  (Test  2  in  Fig.  7,  see 
also  the  data  reported  in  Table  1)  by  the  R(I)  expression  above 
described. 

Due  to  the  choice  of  the  mathematical  relation  to  express  the 
change  of  resistance  with  current,  the  computation  is  carried  out 
with  the  values  of  the  total  cell  resistance  and  current  intensity. 
Then  the  resistance  referred  to  unit  of  geometric  surface  of  MEA 
is  obtained.  It  must  be  observed  that  the  term  R(I)  can  include 
non-ohmic  effects  due  to  the  porous  3D  electrodes.  These  effects 
can  be  related  to  two  different  conditions:  (i)  variable  distribu¬ 
tion  of  ionic  and  electronic  components  of  the  current  within 
the  electrodes  and  (ii)  complexity  of  the  mass  transfer  process 
in  the  porous  mass,  which  may  produce  a  slight  bending  in  the 
middle  part  of  the  polarization  curve,  as  proved  by  many  mathe¬ 
matical  models  accounting  for  diffusion  in  the  porous  mass  [35]. 
However,  these  contributions  are  very  small  for  the  thin  catalytic 
layers  usually  employed  in  PEMFC. 

The  computation  of  the  total  fitting  equation  is  carried  out 
according  to  the  previous  scheme.  To  illustrate  the  interesting 
results  thus  obtained,  the  steps  described  in  the  previous  sections 
can  be  summarized  as  follows: 

1 .  The  linear  fitting  of  the  first  part  of  the  V/log  (|z|)  plot  provides 
the  values  of  the  exchange  current  density  z'o  and  cathodic 
transfer  coefficient  acat  for  oxygen  reduction.  Then  the  acti¬ 
vation  overpotential  for  ORR  is  computed  using  Eq.  (2), 
which  is  numerically  inverted  to  obtain  t]act,O2(0- 

2.  The  open  circuit  voltage  is  evaluated  by  computing  the  mixed 
potential  on  the  basis  an  arbitrary  oxidation  reaction,  Eq.  (7). 
This  allows  the  determination  of  an  estimate  of  the  hydrogen 
crossover  current  density. 

3.  The  value  of  the  exchange  current  density  for  HOR  is 
obtained  and  the  corresponding  activation  overpotential  is 
computed  by  Eq.  (2),  which  is  inverted  to  provide  rjac t,H2(0- 

4.  The  fitting  for  the  ohmic  resistance  value  is  carried  out 
employing  the  relationship  R  (7)  =Ro  (1  +  at1). 

5.  Finally  the  concentration  overpotentials  are  computed,  based 
on  the  known  values  of  the  stoichiometric  coefficients  for 
hydrogen  and  oxygen. 

Having  computed  all  necessary  terms,  the  V(i)  curve  is 
obtained  according  to  Eq.  (1). 

Fig.  5,  which  refers  to  the  polarization  data  of  Test  3  in 
Fig.  7,  illustrates  one  example  of  such  analysis.  Together  with 
the  experimental  data  and  the  total  fitting  equation  (curve 
e),  also  single  terms  are  plotted  in  sequence  in  Fig.  5  as: 
curve  a:  Erev  -  t)act,02;  curve  b:  Erev  -  (f?act,02  +  '7act,H2); 


and  02  (100%  RH)  at  atmospheric  pressure,  H2  stoichiometry  flow  1.5,  02 
stoichiometry  flow  2.  The  letters  from  a  to  d  refer  to  the  various  computed 
polarizations  (activation  overpotential  for  oxygen  reduction,  for  hydrogen  oxi¬ 
dation.  concentration  overpotential  and  ohmic  drop).  Curve  e  is  the  total  fitting 
equation. 

curve  c:  Erev  -  (t?act,o2  +  7act,H2  +  Vconc,0%  +  7conc,H2);  curve 
d:  Erev  -  (t?act,02  +  7act,H2  +  7conc,02  +  7conc,H2  +  Ro  0;  Curve 
e :  Erev  —  [t?act,02  +  7act,H2  +  7conc,02  +  7conc,H2  +  R(0  X  $]. 

In  order  to  compare  the  catalytic  behavior  of  different  MEAs, 
on  the  basis  of  the  analysis  of  the  polarization  curves,  a  parame¬ 
ter  suitable  for  the  evaluation  of  the  catalyst  activity  for  ORR  has 
been  defined  as  the  current  density  at  the  overpotential  |  rj  |  =0.3  V 
(ijpO.3)-  Similar  definitions  have  been  proposed  in  the  litera¬ 
ture  [36],  the  present  one  does  not  refer  to  a  cell  voltage  value 
(typically  0.9  V),  because  different  temperature  and  gas  pres¬ 
sure  conditions  may  change  the  reversible  cell  voltage,  so  the 
choice  of  an  overvoltage  value  (Vrev  —  V)  appears  more  appro¬ 
priate.  The  proposed  parameter  is  obtained  on  the  computed 
curve  Erev  —  t)act,o2,  where  the  activation  overpotential  for  oxy¬ 
gen  reduction  is  only  considered.  In  this  way  the  estimate  of 
the  catalytic  activity  does  not  depend  on  the  other  dissipation 
terms,  particularly  the  ohmic  resistance,  which  may  markedly 
affect  the  MEA  performance.  For  Test  2  described  in  Fig.  7, 
we  found  a  value  i,,=o.3  =  19.2  mA  cm-2,  as  illustrated  in  Fig.  6. 
Having  examined  numerous  tests  with  laboratory  prepared  and 
commercial  MEAs,  we  may  affirm  that  this  parameter  can  be 
very  effective  in  comparing  the  ORR  catalytic  activity. 

The  analysis  described  in  this  paper  has  been  performed  on 
various  tests  corresponding  to  different  performance  of  labora¬ 
tory  prepared  MEAs  and  commercial  MEAs.  As  an  example 
the  results  obtained  by  analyzing  three  different  tests,  chosen 
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Fig.  6.  Determination  of  the  parameter  (i,=o.3)  assumed  to  quantify  the  catalytic 
activity  for  ORR.  The  data  refer  to  test  2  in  Fig.  7. 


Fig.  7.  Plot  of  three  polarization  curves  obtained  on  different  MEAs,  analysed  in 
the  paper.  Tests  1  and  2  refer  to  laboratory  prepared  MEAs  of  5  cm2 ,  with  a  Pt  load 
of  0.5  mg  cm“2,  humidified  H2  and  02  (Test  1 : 60%  RH  and  Test  2: 100%  RH)  at 
atmospheric  pressure,  stoichiometry  flow  2 .  Test  3  refer  to  a  commercial  MEA  by 
ElectroChem  Inc.,  5  cm2, 1  mg  Pt  cm-2,  T=  70  °C,  humidifiedH2  and  02  (100% 
RH)  at  atmospheric  pressure,  H2  stoichiometry  flow  1.5,  02  stoichiometry  flow 
2. 

because  of  markedly  different  performances  (Fig.  7)  are  shown 
in  Table  1.  All  three  tests  were  carried  out  with  MEAs  of  5  cm2. 
Tests  1  and  2  refer  to  laboratory  prepared  MEAs  (as  indicated 
in  the  experimental  section)  with  a  Pt  load  of  0.5  mg  cm-2,  the 
resistance  value  for  Test  1  is  very  high,  compared  to  the  one  of 
Test  2,  because  this  test  was  performed  under  low  humidifica¬ 
tion  conditions  (60%  RH).  Test  3  refer  to  a  commercial  MEA 
by  ElectroChem  Inc.,  with  1  mg  cm-2  of  Pt. 

4.  Conclusion 

A  simplified  model  is  proposed  to  analyse  polarization  curves 
of  single  PEMFC,  by  computing  the  different  overpotential  com¬ 
ponents  which  provide  the  best  fitting  of  the  experimental  data. 

The  open  circuit  potential  for  the  oxygen  electrode,  which 
is  known  to  be  markedly  lower  than  the  reversible  potential,  is 


interpreted  in  terms  of  mixed  potential,  assuming  an  arbitrary 
parasitic  oxidation  reaction.  This  allows  a  very  good  fitting  of  the 
polarization  data  at  extremely  low  current  densities.  In  addition, 
from  the  value  of  the  current  density  of  the  parasitic  reaction 
at  the  open  circuit  potential,  a  good  estimate  of  the  hydrogen 
crossover  can  be  obtained.  This  is  confirmed  by  literature  data. 

A  parameter  for  the  evaluation  of  the  electrocatalytic  perfor¬ 
mance  of  the  ORR  is  introduced.  It  is  defined  as  the  value  of  the 
current  density  for  a  cathode  overpotential  of  0.3  V  (i,,=o.3)  com¬ 
puted  on  the  theoretical  curve  for  the  activation  overpotential  of 
the  ORR. 

The  possible  change  of  the  ohmic  resistance  with  current 
intensity  is  accounted  for  and  discussed. 

The  results  of  three  polarization  tests  performed  on  lab¬ 
oratory  prepared  and  commercial  MEAs  exhibiting  different 
performances  have  been  analysed  with  the  proposed  method, 
which  appears  to  be  an  effective  tool  for  the  evaluation  of  MEA 
performance. 
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